Mon. Not. R. Astron. Soc. 000, [TH6] (2011) Printed 3 November 2011 



(MN style file v2.2) 



Identifying Non- Resonant Kepler Planetary Systems 



Dimitri Veras^*, Eric B. Ford^ 

^Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 OHA, UK 

^Astronomy Department, University of Florida, 211 Bryant Space Sciences Center, Gainesville, Florida 32111, USA 



Accepted 2011 November 01. Received 2011 October 26; in original form 2011 October 09 



o 

(N 

> 
O 



6 



ABSTRACT 

The Kepler mission has discovered a plethora of multiple transiting planet can- 
didate exosystems, many of which feature putative pairs of planets near mean 
motion resonance commensurabilities. Identifying potentially resonant systems 
could help guide future observations and enhance our understanding of planetary 
formation scenarios. We develop and apply an algebraic method to determine 
which Kepler 2-planet systems cannot be in a lst-4th order resonance, given 
the current, publicly available data. This method identifies when any potentially 
resonant angle of a system must circulate. We identify and list 70 near-resonant 
systems which cannot actually reside in resonance, assuming a widely-used for- 
mulation for deriving planetary masses from their observed radii and that these 
systems do not contain unseen bodies that affect the interactions of the observed 
planets. This work strengthens the argument that a high fraction of exoplanetary 
systems may be near resonance but not actually in resonance. 
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1 INTRODUCTION 

The Kepler mission has identified more planetary candi- 
dates than the total number of extrasolar planet s that 
were known before its launch (|Borucki et al.|[201ll ). This 
windfall of mission results enables investigators to per- 
form better statistical analyses to more accurately esti- 
mate properties of the exoplanet population. A dynam- 
ical question of particular interest is, "what is the fre- 
quency of resonant extrasolar systems?" 

The absence of resonance may be just as impor- 
tant, especially in systems with tightly p acked but sta- 
ble configurations (|Lissauer et al.l [201181 ) and/or which 
might feature transit timi ng variations th a t span sev- 
eral orders of magnitude (|Ford et al.l 1201 ll : IVeras et al.l 
I2OIII I. Resonances may be a signature of a particular 
mode of planetary formation; they are likely to be in- 
dicative of convergent migration in nascent protoplane- 
tary disks ( Thommes fc Lissauer" 2003l : iKlev et al.ll2004l : 
|Papaloizou & Szuszkicwicz 2005). Their absence might 
be indicative of a dynamical histo ry dominated by grav i- 
tational planet-planet scattering (|Ravmond et al.ll2008l ). 
Alternatively, resonant statistics might best constrain 
how disk evolution and gravitational scattering interface 
(iMatsumura et al.|[2010: M oeckcl fc Armitage 20ll]). 

Lissauer et al.l ( 2011lJ ) provide a comprehensive ac- 
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counting of all planetary period ratios in /('ep/er systems, 
and present distributions and statistics linking the pe- 
riods to potentially resonant behavior. Our focus here 
simply is to identify which systems cannot be resonant. 
A variety of mean motion resonances are known to exist 
in the Solar System and extrasolar systems. However, for 
most of these cases, the masses, eccentricities and lon- 
gitude of pericenters have been measured directly. Con- 
trastingly, for Kepler systems, these parameters are con- 
strained poorly, if at all. Hence, analyzing resonant Ke- 
pler systems poses a unique challenge. 

In order to address this situation, we present an anal- 
ysis that i) treats the entire potential eccentricity phase 
space in most cases, ii) considers a relevant limiting case 
of the unknown orbital angles, and iii) pinpoints the man- 
ner in which planetary masses can affect the possibility of 
resonance. The analysis is also entirely algebraic, allow- 
ing us to avoid numerical integrations and hence investi- 
gate an ensemble of systems. The results produce defini- 
tive claims about which systems cannot be in resonance, 
given the current data. We consider the ensemble of Ke- 
pler two-planet systems which may be near lst-4th order 
eccentricity-based resonances, under the assumptions of 
coplanarity, non-crossing orbits, and the nonexistence of 
additional, as-yet-unidentified planets that would affect 
the interactions of the observed planets. In Section 2, we 
explain our analytical model, and then apply it in Section 
3. We present our list of non-resonant Kepler systems in 
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Table 1, and briefly discuss and summarize our results in dzu^ _ ^~3^ ^ dRw ^.^^ 
Section 4. di i9e„ ' 

where 



2 ANALYTIC MODEL 
2.1 Proving Circulation 

The behavior of linear combinations of orbital angles 
from each of the two planets determines if a system is 
in resonance or not. "Libration" is a term often used 
to denote oscillation of a resonant angle, and "circula- 
tion" often refers to the absence of oscillation. Resonant 
systems must have a librating angle; for mean motion 
resonances, the focus of this study, the librating angle in- 
corporates the positions of both planets. This seemingly 
simple criterion, however, belies complex behavior often 
seen in real systems^ 

For example, iBrasser et al.l l|2004l 'l illustrate how 
the resonant angle of a Neptunian Trojan can switch 
irregularly between libration and circulation, while 
[Farmer fc GoldreichI (|2006l ) illustrate how long periods of 
circulation can be punctuated by short periods of libra- 
tion. Characterizing libration versus circulation in extra- 
solar systems - often with two massive extrasolar plan- 
ets - sometimes nec essitates careful statistical measures 
l| Veras fc FordlbOld ;). 

The gravitational potential between two coplanar 
bodies orbiting a star can be described by two "disturb- 
ing functions," denoted by R^, with w = 1,2, which are 
infinite linear sums of cosine terms with arguments of the 
form 



(1) 



where A represents mean longitude, represents longi- 
tude of pericenter, the "j" values represent integer con- 
stants, and the subscript "1" refers to the outer body 
while the subscript "2" refers to the inner body. We 
choose this subscri pting convention so our equ atio ns con- 
form w ith those in [Murray fc DermottI l|l999f ) and I Verai 
The mean longitude is directly proportional to 
mean longitude at epoch, denoted by e^, such that 

Au, = I^iu + A/,„ = TVw + + Mwa + Jig nn,{t')dt' = 

Ell) + // y.w^^'' ajuit'Y"^^^^ dt' , where Mm denotes mean 
anomaly, -n-u, denotes argument of pericenter, and n^, de- 
notes mean motion. The mean motion is related to its 
semimajor axis, a^, and mass through Kepler's third law 
by n^a^ — where fiw ~ G (mo -I- rriiu), with mu, rep- 
resenting the planet's mass, mo the central body's mass, 
and G the universal gravitational constant. 

The time derivative of Equation yields: 



^{t) = ji^i'laj(t)+jiii{t)+j2^ilaj(t)+j2e2{t) 



(2) 



Lagrange's Plane tary Equations (e.g., 
iMurrav fc DermottI Il999l . pp. 251-252) relevant to 
Eq. ^ are: 



ddw ^ . QRw , ~ A A dRw /ij-, 

at daw ocw 



(5) 
(6) 

Aw,3 = l^{l~ely. (7) 

The form o f Rw used dictat es how t o proceed. IVerad 

|2003) expresses lElhs fc Murravl 's (|2000l ) disturbing func- 
tion as 



= 2^„ 



y = l 



.P=i 



(8) 



where Cw'^'' is a function of the masses and semimajor 
axes, and X^^'^^ is a function of the eccentricities and 
inclinations. Both of these auxiliary variables contain the 
detailed functional forms needed to model individual res- 
onances. 

When inserted into Eqs. ©-Q, the disturbing func- 
tion will yield: 



a{") 



E( 



D 



+ E 



(9) 



where <j>^^^ is a potentially resonant angle, D^^'"' is an 
explicit algebraic function of the masses and orbital ele- 
ments, and 



E 



= 5Z i^""V™a„ = . (10) 

w — l 

Hence, a particular angle (^'"' cannot librate if 



(11) 



Equation (|lip represents the criterion for an angle 
to circulate. Note that no integrations are required to 
evaluate the criterion. For systems which cannot be in 
resonance, the maximum value of Yl l-D^^'"' l/l-^'"' I = 
provides an estimate of the proximity to a potential res- 
onance. 

The dependence these variables have on planetary 
mass is important to understand for the implications for 
Kepler systems: /ii or fi2 appears in each of the terms 
m B'"^ and Dtf'"), and are insensitive to the planetary 
masses, as long as these masses are negligible compared 
to the star's mass. Further, each term in D^f is linearly 
proportional to either mi or m2. Therefore, if mi and m2 
are both scaled by the same factor of their radii in a com- 
mon mass-radius relationship, then the relative strength 
of the terms in D^'^'"' will vary by a factor proportional to 
the planetary radii ratio. The signs of these terms depend 
on the resonance being studied. Hence, for a particular 
resonance, one may determine which terms are additive, 
and obtain an explicit dependence on planetary mass. 

2.2 Hill Stability 

Unless two bodies are resonantly locked in a configura- 
tion that allows for crossing orbits (such as Neptune and 
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Figure 1. Hill stability eccentricity portrait. Each curve (for 
the following planet/star mass ratios: solid lines - IO^^ttiq; 
dotted lines - lO~*m0; dashed lines - lO~^m0; dot-dashed 
lines - lO~^m0) bounds eccentricities at which two planets 
may be Hill stable, for semimajor axes ratios corresponding 
to the following commensurabilities (in order of dashed lines 
going outward from the origin): 7:6-magenta, 3:2-green, 5:3- 
red, 2:l-blue, and 3:l-salmon. Stable systems lie below the 
curves, and observed Kepler systems are assumed to be stable. 
Note that the curves are not symmetrical about the origin. 
The solid salmon curve indicates that nearly all planets whose 
semimajor axis ratio is within a factor of ~ 2 are subject to 
stability constraints for ei,e2 < 0.3. 



Pluto), the bodies will undergo dynamical instability if 
their eccentricities are too great and their semimajor axis 
difference is too small. Planets who se orbits never cross 
are said to be Hill stable, and obey (|Gladmanlll993l ): 



1 + rii +rj2 



+ riirj2 + ri2) 

V2 



<2 \ 772] 



e^) + 7?iVl 



> 



33r?ir;2 »?i»?2 (11??2 + Tryi) 



(ryi +772)-^ 



3(771+772) 



(12) 



where 771 = 7ni/mo, 772 = 7712 /mo, and a = a2/a\. 
In principle, systems which do not satisfy Eq. (|12|) 
may be stable, but generally this is true only for res- 
onant systems. Thus, we focus our efforts only on 
those systems which are provably stable. Fig. [1] dis- 
plays level curves of Eq. (|12p for different values of 771 
and ?72 (lO~^m0-solid, lO~'*7n0-dotted, lO~^m0-dashed, 
lO~®m0-dashed-dot) and at different commensurabilities 
of interest (7:6-magenta, 3:2-green, 5:3-red, 2:l-blue, 3:1- 
salmon), which each correspond to the appropriate value 
of a. The range of eccentricities plotted and considered 
in this study is well within the a bsolute convergence 
limits of this disturbing func tion (|Ferraz-MeIld 1 19941 : 
ISidlichovskv fc Nesvornv[il994 ) and corresponds r oughly 
to a regime where a fourth-order treatment (as in I Veraa 
|2007| ) is accurate to within « 0.3"* < 1%. 



3 APPLICATION TO KEPLER SYSTEMS 
3.1 Method 

We use the publicly available Kepler databas^B, which 
provides directly measured values for planetary periods 
and estimated values of the stellar mass. With these, we 
estimate the planetary semimajor axes by assuming val- 
ues for the planetary mas s es base d on a mass-radius re- 
lationship. iLissauer et al.l l|2011bl ) estimate these values 
using: 



(13) 



where b — 2.06. A similar formula holds for the other 
planet in the system. 

We assume Eq. (|13|) to be true, and then determine 
the "nearness" of all two-planet systems to a mean mo- 
tion commensurability through 

(1 - x/100) {hlhf' < «i/«2 < (1 + x/WO) (ji/j2)'''' , (14) 

where x effectively measures the percent offset from res- 
onance in semimajor axis space. Because we consider all 
potential resonances up to fourth-order, the same plane- 
tary system may be close to several resonances. We refer 
to a "near-resonant instance" as a case where Eq. H14|) 
is satisfied. We plot the cumulative frequency of near- 
resonant instances as a function of x in Fig. (2] Not in- 
cluded in the figure are (ji , 72) multiplicities: higher-order 
harmonic angles which can be expressed as a linear com- 
bination of the angles with the lowest-order relatively 
prime values of ji and j2, although these higher-order 
terms are i ncluded in t he computation of D^^'"' (see, e.g. 
Table 4 of IVera s 2007). The figure demonstrates the po- 
tentially wide variety of resonant configurations which 
may be possible depending on one's definition of "near- 
resonant". We choose to be conservative and include all 
"near-resonant instances" on the plots in our analysis. 

Having identified a potential resonance with a given 
ji and j2, and having adopted mi, m2, oi and 02, we then 
sample the entire eccentricity phase space for which the 
planets are Hill stable. In the few cases where Hill stable 
systems admit ei > 0.3 and/or €2 > 0.3, we limit the ec- 
centricities to these values for accuracy and convergence 
considerations, as described earlier. At each of 31 evenly- 
spaced values of 62 from to 0.3, we compute the maxi- 
mum Hill stable value of ei, and then sample 31 evenly- 
spaced values of ei from to that Hill maximum value. 
Finally, we compute /3„ = 1^*^'"' l/l-^'"^ | for every 
possible disturbing function argument for ji,|j2| > 
which satisfy the dAlembert relation^. If no value of /3u 
exceeds unity, then we flag the system as non-resonant 
and record the highest value of Pu as a proxy for closest 
proximity to resonance. 



^ at |http: / /archive.stsci. edu/kepler/planet_candidates.html| 
^ The secular arguments ji = J2 = up to 4th-order are 
included in the computation of D^"'"' even though we test 
only for circulation of angles which can lead to moan motion 
resonance. 
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Figure 2. Number of near-resonant instances in 2-planet Ke- 
pler systems. Plotted is the cumulative number of triplets 
(KOI,ji,j)2) vs X, where x is defined in Eq. II14I I and KOI 
refers to the Kepler Object of Interest, or candidate exoplanet 
host star. The upper panel plots all instances for lst-4th order 
resonances, and the other panels plot selected resonances of 
interest. 



3.2 Results 

Of the 116 Kepler systems with two transiting candi- 
dates for which public data is available, 94 have planets 
with periods similar enough to each other to be consid- 
ered "near-resonant" for at least one lst-4th order res- 
onance with X ^ 10. These 94 systems registered a to- 



tal of 465 near-resonant instances for ji < 20. Of those 
465, 313 cannot be resonant. Only 4 systems (KOI #s: 
89, 523, 657, 738) have planets which may be resonant 
for at least one instance and cannot be resonant for at 
least one other, thereby restricting the system's potential 
membership in some resonances. For 20 systems (KOI 
#s: 82, 111, 115, 124, 222, 271, 314, 341, 543, 749, 787, 
841, 870, 877, 945, 1102, 1198, 1221, 1241, 1589), mean 
motion resonance could not be ruled out for any of the 
near-resonant instances sampled. Quantifying the likeli- 
hood of resonance in these systems would require future 
detailed individual analysis. 

We claim that the remaining 70 systems cannot be 
in resonance, despite their close proximity to a commen- 
surability. We list these systems in Table 1, along with 
the number of their near-resonant instances, their maxi- 
mum value of P and the mean motion resonance this value 
corresponds to. The larger the value of /?, the closer the 
system is to having parameters which could admit reso- 
nance. In most cases, /3 < 0.5, meaning that the systems 
are well outside of resonance. Notable exceptions are KOI 
244, KOI 431, KOI 508, KOI 775, and KOI 1396, three of 
which are close to the strong 2:1 commensurability. KOI 
1151 is so close to so many commensurabilities because 
for that system, 1/a — 1.26, meaning that the planets 
are tightly packed and on the edge of stability. 

We don't expect variations in the masses of the Ke- 
pler planets to greatly affect the composition of Table 
1, assuming that the value of b from Eq. p3[) does not 
vary by more than a few tenths from 2.06. The value of 
P for a particular system can hint at the potential impli- 
cations of mass variation. In particular, values of /3 close 
to unity indicate that the system is on the border of po- 
tentially resonant behavior. For example, for KOI 1396 
{P — 0.89), if we set 6 = 1.7, then the resulting plane- 
tary masses could allow the system to harbor a (weak) 
9:5 resonance. 



4 CONCLUSION AND IMPLICATIONS 

We have identified 70 Kepler 2-planet near- 
commensurate systems which cannot be in an 
eccentricity-based mean motion resonance of up to 
4th-order. These systems, may, in principle, achieve reso- 
nance with crossing orbits or high (e > 0.3) eccentricities 
that could remain stable. The criterion of Eq. (jlip is 
generally applicable to any 3-body system suspected 
of harboring resonant behavior. We caution that these 
results could be affected by the presence of additional 
planets which have yet to be detected. 

Systems which are provably non-resonant may be 
incorporated in formation studies and detailed analy- 
ses of Kepler data. Kepler multi-planet systems may 
be preferentially clust ered around particular resonances 
l|Lissauer et al.l2011bl ). One possible explanation may be 
that convergent migration locks planets in a mean motion 
resonance which later gets broken by some additional per- 
turbation. We find that Kepler multi-planet systems pref- 
erentially cluster around commensurabilities where |j2| 
is low and rarely do so when \j2\ is high. In particular, 
the number of Kepler systems near the 2:1, 3:2, 3:1 and 
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Table 1. Non-Resonant 2-Planet Kepler Systems. 



KOI 


Number of 


Closest 


Maximum 


Number 


Near-Res Instances 




13 


116 


1 


3;1 


0.078 


123 


1 


3:1 


0.018 


150 


1 


3:1 


0.027 


153 


4 


9:5 


0.28 


209 


2 


3:1 


0.22 


220 


6 


9:5 


0.039 


232 


3 


7:3 


0.016 


244 


3 


2:1 


0.67 


270 


3 


3:1 


0.0043 


279 


4 


9:5 


0.17 


282 


1 


3:1 


0.016 


284 


2 


3:1 


0.033 


291 


1 


4:1 


0.010 


313 


3 


7:3 


0.043 


339 


1 


3:1 


0.0043 


343 


2 


7:3 


0.12 


386 


2 


5:2 


0.16 


401 


1 


5:1 


0.031 


416 


1 


5:1 


0.024 


431 


2 


5:2 


0.60 


440 


1 


3:1 


0.040 


446 


6 


9:5 


0.18 


448 


2 


4:1 


0.038 


456 


1 


3:1 


0.048 


459 


2 


3:1 


0.042 


474 


3 


5:2 


0.020 


475 


4 


9:5 


0.23 


497 


1 


3:1 


0.31 


508 


2 


2:1 


0.84 


509 


2 


3:1 


0.042 


510 


3 


7:3 


0.035 


518 


1 


3:1 


0.028 


534 


2 


7:3 


0.37 


551 


3 


2:1 


0.25 


573 


1 


3:1 


0.11 


584 


2 


2:1 


0.083 


590 


2 


4:1 


0.0049 


612 


3 


7:3 


0.21 


638 


2 


3:1 


0.12 


645 


2 


3:1 


0.031 


658 


6 


9:5 


0.10 


672 


3 


5:2 


0.12 


676 


1 


3:1 


0.12 


691 


5 


9:5 


0.11 


693 


5 


9:5 


0.17 


700 


1 


3:1 


0.027 


708 


3 


7:3 


0.036 


736 


2 


3:1 


0.063 


752 


1 


5:1 


0.0049 


775 


2 


2:1 


0.78 


800 


3 


3:1 


0.023 


837 


3 


9:5 


0.035 


842 


2 


3:1 


0.090 


853 


6 


9:5 


0.19 


869 


1 


5:1 


0.051 


896 


3 


5:2 


0.14 


954 


2 


5:1 


0.0059 


1015 


2 


7:3 


0.085 


1060 


2 


5:2 


0.028 


1113 


1 


3:1 


0.030 


1151 


15 


9:7 


0.044 


1163 


2 


3:1 


0.015 



Table 2. ...Table 1 continued 



KOI 


Number of 


Closest 


Maximum 


Number 


Near-Res Instances 


jl-\j2\ 


13 


1203 


3 


7:3 


0.059 


1215 


4 


9:5 


0.084 


1278 


1 


4:1 


0.0081 


1301 


1 


3:1 


0.13 


1307 


3 


7:3 


0.049 


1360 


2 


5:2 


0.27 


1364 


1 


3:1 


0.19 


1396 


6 


9:5 


0.89 



5:3 commensurabilities is higher than what would be ex- 
pected from a random distribution of Kepler planet can- 
didate semimajor axes. Provably non-resonant planets 
may also complement transit timing variation statistics, 
as these variations take on distinctly different characteris- 
tics f or near-resonant and resonant systems (|Veras et al.l 

[mH). 

The analysis in this work cannot be performed with 
Kepler systems that contain more than 2 planet candi- 
dates because i) more disturbing functions must be incor- 
porated, and hence the criteria for circulation becomes 
decidedly more complex, and ii) analyt ical formulae for 
Hill st ability no longer hold. Figure 29 of lChatteriee et al] 
(|2008l ) demonstrate that widely-separated pairs of plan- 
ets in 3-planet systems whose orbits would be Ifill stable 
in the 2-planet-only case may eventually become unsta- 
ble. Even if a multi-planet system was assumed to be sta- 
ble over a specified period of time and more disturbing 
functions were introduced, the resulting expansion of the 
phase space might render the computational cost of a sim- 
ilar algebraic analysis prohibitive compared to numerical 
integrations. However, the investigation of the nonplanar 
2-planet case with inclination-based resonances might be 
a fruitful avenue to explore, especially because multiple 
transits detected by Kepler may constrain the planets' 
mutual inclination, albeit weakly. 
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